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We present a detailed microscopic study of quasi-ballistic transport in deep submicron semiconductor 
channels. In particular, we study the crossover between the diffusive and ballistic regimes of transport and 
identify signatures in the electrostatic response, electron density, and nonequilibrium electron distribution 
function that are due to ballistic and diffusive transport, respectively. Our theoretical and computational 
approach is based on the Boltzmann transport equation for a nondegenerate electron system, together with 
the Poisson equation, from which we obtain the nonequilibrium electron distribution in a self-consistent way. 
We find that the electron distribution is significantly different from a near-equilibrium, shifted Maxwell- 
Boltzmann distribution function, and that it displays a large broadening, as well as pronounced features, 
peaks and shoulders in the high- velocity tail of the distribution. These features are signatures of the nonequi- 
librium and quasi-ballistic nature of the electron transport. An analysis of the spatial density distributions 
of electrons in the semiconductor channel region shows that the relative scattering efficiency varies within 
the channel region. Our study of the crossover between the diffusive and ballistic regimes of transport, which 
we investigate by downscaling the size of the semiconductor channel, shows that a crossover occurs around 
a channel length specific to the material parameters, and is clearly accompanied by a redistribution of the 
voltage drop between the contacts and the channel region. For these transition channel lengths, we also 
observe clear signatures of the diffusive-to-ballistic crossover in the spatial electron density distribution that 
manifest as exponential-like and linear spatial dependences, indicative of diffusive and ballistic transport, 
respectively. 
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I. INTRODUCTION 

Advances in fabrication and processing of semiconduc- 
tor materials on the nanoscale are presently enabling the 
realization of deep submicron and nanoscale semiconduc- 
tor structures and devices. At these ultrasmall length 
scales, finite size effects, hot-electrons, ballistic transport 
and microscopic doping variations are important in de- 
termining transport properties. In addition, these sys- 
tems can be very easily brought out-of-equilibrium, as 
very high fields are established even for low applied bias 
voltages. 

The physics of electron transport in ultrasmall semi- 
conductor systems is thus very rich as well as com- 
plex. Drift-diffusion approaches as well as approxima- 
tions based on a drifted-Maxwell-Boltzmann distribution 
break down and an adequate treatment of the nonequi- 
librium distribution function is needed. In particular, 
the high-energy tail of the distribution can have com- 
plex features that are determined by the scattering pro- 
cesses and the electrostatics within the systems as shown 
experimentally in Raman spectroscopy measurements of 
the electron distribution in submicron, inhomogeneous 
III-V systemsi These experiments have shown far-from- 
cquilibrium features such as strongly broadened velocity 
distributions and interesting high- velocity tail structured 

Theoretically, there have been a few recent works 
that have shown that high-velocity structure in the 
distribution function can be attributed to ballistic 
effects i 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ^ 1 In a recent paper by us^ we dis- 
cussed the spatial dependence of the distribution func- 
tion in deep submicron GaAs channels and showed how 



quasi-ballistic signatures manifest as a broadening and 
high-velocity tail structure, as seen in the experiments. 
We also analyzed microscopic features of scattering in the 
channel region in the presence of energy-dependent scat- 
tering due to polar-optical phonons. Interesting similar 
analyses with or without the inclusion of finite scatter- 
ing have also been previously reported within the same 
theoretical frameworks 

The scattering problem in submicron and nanoscale 
channels is complicated since scattering and the electro- 
static response of the system constitute a strongly cou- 
pled physical problem that creates a self-consistent feed- 
back and an interesting charge redistribution. Since pure 
ballistic transport is not attainable at room temperature 
in the presence of finite scattering 3 ^ and thus, the trans- 
port in submicron systems is quasi-ballistic to a varying 
degree, it is interesting to ask: What are the signatures of 
nonequilibrium and quasi-ballistic transport in quantities 
other than the electron distribution, such as electron den- 
sity, potential energy and how can one separate the ballis- 
tic contribution to the problem and study its influence on 
the total quantities? More importantly, what constitutes 
the crossover between diffusive and quasi-ballistic trans- 
port and how does it manifest on a microscopic level? 

In this paper we report on a detailed study of 
the crossover between the diffusive and quasi-ballistic 
regimes of transport in (sub)micrometer scale semicon- 
ductor channels, in particular investigating how this 
crossover manifests in the electrostatic response of the 
system, the electron redistribution within the channel 
region, the contribution of ballistic electrons to the to- 
tal electron density, and the nonequilibrium electron dis- 
tribution function. Our theoretical and computational 
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framework is based on a self-consistent solution of the 
Boltzmann and Poisson equations, 12 which enables us to 
calculate nonequilibrium electron distributions and cor- 
responding moments, as well as the accurate electrostatic 
reponse. 

We identify that, for the GaAs system that we study, 
with a given doping profile and at room-temperature, a 
crossover to quasi-ballistic transport occurs for channel 
lengths smaller than w 1.5 ^m. In the electrostatics, the 
crossover manifests as as redistribution of the potential 
drop from the channel to the contact regions, for decreas- 
ing channel lengths. In the electron density, which we 
study by separating components with different velocity, 
we observe interesting spatial dependences reflecting the 
diffusive-to-ballistic crossover. Finally, on a microscopic 
level, we examine the nonequilibrium electron distribu- 
tion function and its phase-space dependence around the 
crossover regime. The presented results provide insight 
into the role of scattering, ballistic transport and elec- 
trostatics in the transport properties of deep submicron 
and nanoscale semiconductor systems. 

The paper is organized as follows: First, we give a brief 
introduction to the theoretical model in Section^ refer- 
ring the reader to Ref. [l2| for further details. Next, we 
will present our results and discussion, first presenting 
an introduction to the basic signatures of quasi-ballistic 
transport in Section IIII Al Subsequently, we discuss 
the crossover between the diffusive and quasi-ballistic 
regimes of transport in Section 1111 Bl We conclude and 
summarize our main results in Section llVI 



II. MODEL 

A microscopic study of the kinetics of quasi-ballistic 
transport in ultra-small inhomogeneous semiconductor 
structures requires in principle the solution of the Boltz- 
mann transport equation (BTE). In addition, in order 
to capture the full physical picture, the BTE needs to 
be solved self-consistently with the Poisson equation. 
This is a very difficult task since the BTE is a non- 
linear, integro-differential equation for the semiclassical 
distribution function, /(r,v,i), which in principle has 
seven dimensions. Although the Monte Carlo method 
has been very popular and successful for the solution of 
the BTE in semiconductor device simulationji&ii sev- 
eral worksiSii^iiii have recently solved the BTE by direct 
methods. This allows noise-free spatial and temporal res- 
olution of the electron distribution function, which in the 
Monte Carlo method may be difficult to obtain due to the 
statistical nature of the approach. 

In our model, the electron transport is described by the 
following one-dimensional BTEli within the relaxation 
time approximation 

eE(x)df(x,v) df(x,v) f(x,v)-f LE (x,v) 



dv 



dx 



where m* is the effective mass in the parabolic band 
approximation and where we have assumed that the 
nonequilibrium electron distribution relaxes to a local 
equilibrium distribution, /l£(.t,i>), which we in the fol- 
lowing assume to describe non-degenerate statistics ac- 
cording to the following normalized Maxwell-Boltzmann 
(MB) distribution 
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where n(x) is the local electron density, and T is the lat- 
tice temperature. The electron density is obtained from 
the full electron distribution according to 
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The BTE, eq. is coupled to the Poisson equation 

through the electric field and electron density. In the fol- 
lowing we assume that our system is unipolar and hence, 
that the Poisson equation is given by 
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where Njj(x) is a given inhomogeneous doping distribu- 
tion profile and e is the dielectric constant. 

Equations 111 II) are coupled through the electron den- 
sity and electric field, and thus, they need to be solved 
self-consistently. We use a numerical approach based on 
finite difference and relaxation methods, the details of 
which can be found in Ref. 0. As boundary conditions, 
we adopt the following scheme: For the potential, the 
values at the system boundaries, £;( e /t) an d x r ^ ight ^, are 
fixed to 4>{ x i) = Vb an d 4>{x r ) = 0, respectively, corre- 
sponding to an externally applied voltage Vb- The elec- 
tron density is allowed to fluctuate around the system 
boundaries subject to the condition of global charge neu- 
trality, which is enforced between each successive itera- 
tion in the self-consistent Poisson-Boltzmann loop. The 
size of the highly-doped contacts is chosen to be large 
enough such that the electron density is constant deep 
inside the contacts. In addition, the size of the contacts 
has to be large enough, such that the electric field deep 
in the contacts is constant, and very low. This allows us 
to use the analytical, linear response solution to the BTE 

f(x l<r ,v) = fLE(xi, r ,v) [1 - vE(xi r )r / kgT] , (5) 

as phase space boundary conditions at xi. r , where we 
use the local equilibrium distribution, fLE(%i,r,v) and 
local electric field, E(x^ r ), obtained from the previous 
numerical solution to the Poisson-Boltzmann iterative 
loop. At the velocity cut-off in phase space, we choose 
f(x, Umax) = /(i,-«max) = !le{x,v), which is reason- 
able since, in the calculations, we assume t> max > 30/csT. 
A more detailed description and discussion of our numer- 
ical method can be found in Ref. [l2j 
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III. RESULTS AND DISCUSSION 

In the following, we will present and discuss our 
main results. First, we will discuss some basic signa- 
tures of quasi-ballistic transport, followed by a study 
of the crossover between the diffusive and quasi-ballistic 
regimes of transport. Our model structure consists of a 
lightly doped region (hereafter called the channel), with 
doping concentration N — 10 19 m -3 , sandwiched be- 
tween two highly doped semiconductor contact regions 
with doping concentration N + = 10 23 m~ 3 . We assume 
that the entire N + — N — N + structure consists of GaAs, 
and use the appropriate parameters in the calculations 
(effective mass to* — 0.067too, where too is the bare 
electron mass, and dielectric constant e = 13.1). Unless 
otherwise indicated, we assume that the lightly doped 
channel region is 200 nm long. In addition, all of our cal- 
culations are performed at T = 300 K, using a scattering 
time of t = 0.25 ps, which corresponds to a realistic 
scattering time due to polar optical phonon scattering^ 
The microscopic effects of an energy-dependent scatter- 
ing mechanism have been previously studied in Ref. 0- 
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A. Basic signatures of quasi-ballistic transport 

The description of transport in submicron, inhomo- 
geneous semiconductor structures such as the ones de- 
scribed in our paper is a complicated process due to i) the 
size of the system, ii) the inhomogeneous field distribu- 
tion, in) the strongly non-equilibrium and quasi-ballistic 
nature of the problem, iv) the self-consistent feedback 
between the scattering processes and electrostatics. In 
the following, we will discuss some general properties of 
our model N + — N — N + structure, which are valid for a 
wide range of ultrasmall, inhomogeneous semiconductor 
structures and devices. Some of these results have been 
previously discussed by us in Ref. U 

In Fig. da) the electron potential energy (solid line) 
and electric field (dotted line) are shown for a N + — 
N — N + structure with dimensions 2.9/0.2/2.9 /jm, for 
an applied bias voltage of V& = —0.3 V. The potential 
energy has a very slow spatial dependence in the large 
contacts resulting in only ss 50% of the voltage drop 
occuring within the contact regions. In contrast, most 
of the potential drops over the lightly doped, deep sub- 
micron channel region, resulting in a dramatic spatial 
variation. Correspondingly, a large field build-up, with 
opposite sign, occurs at the N + /N interfaces, a strongly 
inhomogeneous electric field is observed in the channel 
region, and in contrast, the field in the contacts is very 
small and constant, thus validating the assumption of 
a near-equilibrium solution to the BTE at the system 
boundaries. The observed strong spatial dependence of 
the electron potential energy and electric field in the 
channel region is due to the carrier diffusion induced by 
the large variation in doping concentration. 2 

In addition to the electrostatic picture discussed above, 
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FIG. 1: (Color online) (a) Electron potential energy and 
electric field for a N + — N — N + structure with doping 
10 23 /iu 19 /i0 23 m - 3 for an applied bias vo itage of V b = -0.3 

V, temperature T = 300 K, and with scattering time r = 
2.5 • 10~ 13 s. Vertical dashed lines mark the interfaces be- 
tween the highly doped N + and lightly doped N regions, (b) 
Electron densities, n + , n~ , n, and (n + — n~)/(n + + n~) as 
defined in the text, (c) "Ballistic" electron densities, tig, n B , 
and "ballistic" electron fractions, n^/n, and rig/n as defined 
in the text. 

valuable information can be obtained from the carrier dis- 
tribution within the system for a given applied external 
field and inhomogeneous doping profile. In this context, 
it is interesting to study the separate contributions, n + 
due to electrons with v > (hereafter denoted as v + ), 
and n~ due to electrons with v < (hereafter denoted 
as v~), to the total electron density 

n{x) — n + {x) + n~ (x) — / f(x,v)dv+ / f(x,v)dv. 

J 17+ J V~ 

(6) 
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In addition we also introduce and study the "ballistic" 
electron density, 

n B (x) = / f(x,v)dv , (7) 

J \v \>v p 

which corresponds to the density of electrons with ener- 
gies higher than the top of the potential barrier between 
the source contact and the channel [occurring at x-i in 
Fig. da)], thus corresponding to a "ballistic" electron 
density. In eq. Q, v p = y/2/m*\Ua — U(x)] and Uq 
and U (x) are the barrier top and local electron potential 
energies, obtained from the Poisson equation Q, respec- 
tively. 

In Fig. db) the electron densities n + (dashed line), 
n~~ (dotted line), and n = n + + rT (solid line), i.e., the 
total electron density, are shown for the structure with 
the same parameters as discussed in Fig. d a ) • I n addi- 
tion, we show the spatial dependence of the surplus of 
electrons with velocities v + relative to the total electron 
density, P n = (n + — n~)/{n + + n~) (dashed-dotted line) , 
which can be viewed as a "polarization" of the electrons 
in velocity space. 

Several observations are made: First, the total electron 
density, n, displays an asymmetric spatial dependence 
with respect to the direction of transport, indicating a 
pile-up of charge at the source-channel interface. Second, 
it is seen that the difference between the two densities, 
n + — n~ , has a significant spatial dependence, and that 
the minimum of the electron density components occurs 
at different points along the channel. Furthermore, P n 
attains its maximum around mid-channel. 

One can view the above results within the following in- 
tuitive picture: As electrons are injected from the source 
contact into the channel, they scatter at a rate given by 
1/t. However, the mean free path of the v + and v~ 
electrons is significantly different due to the large built- 
in electric fields around the channel region, giving rise 
to a considerable enhancement of the n + component of 
the density in comparison to the n~ . The latter on the 
contrary is suppressed since, i) electrons injected from 
the drain have to overcome a large potential barrier, ii) 
source-injected electrons scattered from the n + contribu- 
tion and back to the n~ contribution of the density have 
a relatively short mean free path due to the built-in field 
that exerts a force in the opposite direction and thus, 
quickly scatter back to the n + component. The spatial 
dependence of n + , n~ , and (n + — n~)/(n + + n~) indi- 
cates that the efficiency of the scattering is non-uniform 
within the channel. 

It is interesting to analyze the previously defined "bal- 
listic" electron density, ns, and compare with the above 
shown results for the electron densities calculated from 
the entire velocity range. In Fig. d c ) the thick solid 
and dashed lines correspond to the v + , and v~ compo- 
nents (which we in the following denote by ri^ and rig) 
of ns , respectively. In addition, we also show the "ballis- 
tic" electron fractions n~^/n (thin dotted line), and rig/n 
(dashed-dotted line). 



To the left of the source-channel interface at x = —0.1 
/im, the two components, and rig, have a slow spa- 
tial dependence and are only slightly different in value, 
reflecting the small field that is present in the source 
region. Around x — x%, corresponding to the top of 
the potential barrier [see Fig.d a )L the two components 
rapidly split, and peaks at the barrier top, as expected 
by the pile-up of source-electrons at the interface, where 
interestingly, only «80% of the total density consists of 
electrons with positive velocity as seen in n^/n. Since 
the potential barrier seen by the drain-injected electrons 
from the right is far too large compared to ksT to explain 
the large rig component at X2 in terms of thermionically 
injected electrons from the drain, the main contribution 
to that component consists of backscattered source in- 
jected electrons. This is consistent with the interesting 
analyses by Sano&* in which it was shown mathemati- 
cally that, due to the singular nature of the BTE (note, 
that the analyses was made for the BTE in the relaxation 
time approximation), a finite rig around the top of the 
barrier separating the source and the channel region, is 
always present provided there is a finite scattering rate. 
The intuitive physical picture is that, electrons close to 
the conduction band edge can have an infinitely small ve- 
locity, and thus can spend long enough time to experience 
a scattering event. If electrons incident from the source 
experience a sharp potential rise around the barrier, elec- 
trons pile up close to the conduction band edge and are 
successively slowed down, thus increasing the time spent 
around the barrier region, and correspondingly giving rise 
to a significant contribution in rig. The contribution to 
n~ becomes larger the larger the curvature of the poten- 
tial energy around the barrier region. 4 

Beyond the location of the top of the barrier, for x > 
X2, rig drops rapidly within the channel region, followed 
by a change in the slope past the channel-drain interface. 
The positive component, n't, on the other hand, has a 
much slower decay within the channel region, followed 
by a similar change of decay rate at the channel-drain 
interface as seen in rig . We note that rig' is two orders of 
magnitude larger than rig at the channel-drain interface 
showing that, indeed there are very few drain-injected 
electrons at kinetic energies comparable to the potential 
barrier at the source-channel interface. 

A simple exponential fit to rig' and rig in the channel 
region shows that the mean free paths of electrons with 
v > and v < are rj 140 and 50 nm, respectively. 
Hence, source-injected electrons have a mean free path 
comparable to the channel region and, thus, propagate 
quasi-ballistically. However, inspite of the quasi-ballistic 
transport of these electrons, their relative contribution 
to the total electron density in the drain is very small, 
as shown by the thin dashed and thin dashed-dotted 
line representing rig/n and rig/n in Fig. d c )- As ex ~ 
pected, the fractions rig/n and rig/n reach their maxi- 
mum at the barrier top and subsequently decrease mono- 
tonically in the channel. We note, however, that n~gjn 
has an exponential-like decay rate, indicative of diffusive- 
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FIG. 2: (Color online) (a) Normalized electron distribution, 
f(xi,v)vth/n(xi), for the structure studied in Fig. Q calcu- 
lated at the spatial points xi-xe as marked by the arrows in 
Figs- Da) andUb), for V b = -0.3 V. (b) Contour plot of 
f(x,v)v t h/n{x) around the lightly doped N region. 

like transport, in contrast to n^/n which decays almost 
linearly. 

Full information about the microscopic details of the 
electron transport is contained in the nonequilibrium 
electron distribution function, to which we have access 
through our self-consistent solution of the BTE. In Fig. 
I2a) we show the normalized electron distribution func- 
tion, f{xi,v)vth/n{xi), where v t h is the thermal velocity 
and n(xi) is the electron density at calculated at dif- 
ferent locations: one point in the source contact, four 
points along the channel region, and one point in the 
drain contact [same as arrows in Fig. [Ja)]. The elec- 
tron distribution shows a strong spatial dependence, and 
a far-from-equilibrium form at all but two of the spa- 
tial points depicted in Fig. Ufa), namely the ones in the 
contacts. The thick dashed curve, corresponding to the 
distribution at x\ = —0.15 [atxi, is nearly Maxwellian, 
showing that the electrons at that point are near equi- 
librium. This is in agreement with the fact that X\ is 
situated in the contact region, where the field is very 
small and constant [see Figs, da) an d 13b)] and thus, 
the electrons are distributed according to a shifted MB 



distribution, as given by the linear response solution to 
the BTE, eq. ©. 

The electron distributions at xi — x^ on the other hand 
are highly asymmetric, displaying a strong weight toward 
positive velocities, as well as pronounced structure in the 
high- velocity tail for v + . At X2 — —0.05 /zm, the dot- 
ted curve shows the electron distribution at the top of 
the potential barrier between the source and channel re- 
gions. The distribution is distorted from a Gaussian-like 
MB distribution and shows a suppression of electrons in 
the v < tail of the distribution as expected since, in- 
tuitively, only electrons injected from the source, with 
positive velocities, should be present. The presence of 
a negative-velocity tail in the distribution at the top of 
the potential barrier is again in accordance with our dis- 
cussion above and the findings of SanOf^ showing that 
significant backscattering occurs around the barrier re- 
gion due to the sharp potential barrier and corresponding 
pile-up of electrons at the conduction band edge result- 
ing in an increased backscattering efficiency due to the 
prolongued dwell time in the barrier vicinity. 

Beyond the barrier top, the peak of the distribution 
rapidly shifts toward higher velocities and is significantly 
broadened (x^ = 0). In addition, beyond the channel 
half-point, the low- velocity contribution of the distribu- 
tion function increases, and a second peak emerges, as 
seen in the dash-dotted line, corresponding to X4 = 0.05 
/im. At X5 = 0.09 /im, the low-velocity peak dominates 
and only a shoulder is observed on the high- velocity side 
of the non-equilibrium distribution. Finally, at x§ = 0.15 
/im, i.e., in the drain contact, the distribution function 
again assumes a near-equilibrium form. Note however 
that a high-velocity tail is still observable, and the dis- 
tribution function, although closely resembling a shifted 
MB distribution, is significantly different from the one at 

A contour plot of f(x,v)vth/n(x) around the channel 
region is shown in Fig. H^b) and further illustrates the 
general features. The emergence of a narrow peak and 
its rapid shift toward high velocities for x > X2 are clearly 
seen. The horizontal lines indicate the points x\-xq cor- 
responding to the curves in Fig. EJa) and the arrows in 
the potential energy in Fig. Ofa) ■ 

The observed features are a signature of quasi-ballistic 
transport. The high- velocity peak in the distribution cor- 
responds to ballistic electrons which, initially injected 
over the source-channel barrier, due to the rapid and 
large (compared to ksT) potential drop and correspond- 
ing strong inhomogeneous electric field, are accelerated 
to high (positive) velocities. As discussed before, in the 
process, the mean free path of backscattered electrons is 
strongly reduced by the large electric field, and thus, the 
"ballistic" peak reaches very high velocities, and pene- 
trates deep into the channel region, as seen in Figs.^a) 
and[2fb). However, in addition to the narrow "ballistic" 
peak, a broad distribution of electrons occupy the low 
and negative velocity part of the phase space. This com- 
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ponent of the distribution function becomes more and 
more pronounced beyond the channel mid-point. The 
origin of this broadening at the drain end of the channel 
is due to electrons with v < 0, injected from a near- 
equilibrium distribution in the drain contact which pen- 
etrate the channel and contribute to the negative and low 
velocity components of the distribution function. Within 
the channel, the inhomogeneous electric field causes a 
non-uniform dependence of the mean free paths of the n + 
and n~ components of the electron density, as discussed 
previously, which in turn creates a feedback to the elec- 
trostatics due to the charge (re)distribution. It is thus 
evident from Figs.|2Ja) and^b) that the microscopic de- 
tails of the influence of scattering in a submicron channel 
can be obtained only from the full non-equilibrium dis- 
tribution function, which can only be obtained by the 
self-consistent solution of the BTE. 



B. Crossover between quasi-ballistic and diffusive 
transport 

In the following, we will discuss the crossover between 
the quasi-ballistic and diffusive regimes of transport in 
terms of numerical calculations and a study of the mi- 
croscopic signatures discussed above. We will focus on 
the quantities discussed in the previous section and in- 
vestigate their dependence on an externally applied bias 
voltage, as well as the dependence on the channel length. 
We have performed calculations for different structures 
with channel lengths ranging between 0.2 to 4 /im. In 
order to exclude the influence of finite contact size on 
the results, we have performed and compared studies us- 
ing two definitions of the channel length variations. In 
one case, which we simply call the constant contact (CC) 
size case, we keep the size of the contacts fixed, while 
varying the channel length, thus effectively changing the 
total sample size. The second case, which we simply call 
the constant sample (CS) size case, we assume that the 
total sample size is constant, such that the variation of 
the channel length effectively reduces the contact size. In 
both cases, as will be seen in the presented results, the 
contact sizes are found to be large enough not to influence 
the main features of our study. 



1. Electrostatics 

In Figs. |2| we show our calculated results of the elec- 
trostatics for different applied bias voltages and differ- 
ent channel lengths. In Fig. |^a), the potential energy 
profiles (thin lines) and electric field distributions (thick 
lines) are shown for an N + /N/N + structure with dimen- 
sions 2.9/0.2/2.9 /im and the same parameters used in 
the previous results, calculated at the applied bias volt- 
ages Vb — 0, —0.1, —0.2, —0.3 V. Even in the absence of 
an externally applied field, a large potential barrier and 
a corresponding strong and inhomogeneous electric field 
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FIG. 3: (Color online) (a) Potential energy (thin lines) and 
electric field (thick lines) profiles for different external bias 
voltages, for the same parameters and structure studied in 
Figs. and |2] (b) Fraction of the potential drop that 
occurs across the contacts, V CO ntacts/Vb, and the channel, 
Vchanne,l/Vb, as a function of channel length. Inset shows the 
spatial dependence of the potential energy profile for channel 
lengths in the range L c — 0.2 — 4 /im. The different curves 
correspond to the constant contact (CC) and constant sam- 
ple (CS) size approaches to the channel length variation in 
the calculations. 



builds up in the channel region due to the interdiffusion 
of electrons from the highly doped contact regions to the 
lightly doped channel region. As an external bias is ap- 
plied to the sample boundaries at ±3 /im, the majority 
of the potential drops over the channel region. We note 
that, for the considered bias range, Vb = —0.1 to —0.3 V, 
the fraction of the potential drop that takes place over the 
channel region is effectively the same and « 50%. The 
electric field, although varying withing the same range of 
« —25 to w 15 kV/cm has an increasingly stronger spatial 
dependence with increasing bias voltage. 

To investigate the crossover between the diffusive and 
quasi-ballistic regimes of transport in the electrostatics, 
we have studied the electrostatic response for structures 
with different channel lengths, ranging between L c = 0.2 
and 4 /xm. The potential energy profiles for different 
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channel lengths (calculated for the CS case) at Vb = —0.3 
V are shown in the inset of Fig.[3{b). It is evident from 
the results that, although most of the potential drop, 
and hence resistance, occurs over the channel region, the 
fraction of the potential drop associated with the channel 
region strongly depends on the channel length. 

In Fig. OUb) we show the fraction of the potential drop 
that occurs across the channel, V channel I Vb-, and contact, 
Vcontacts/Vb: regions, respectively, for different channel 
lengths. For long channels, virtually all of the applied 
voltage drops across the channel region, irrespectively of 
the magnitude of the applied voltage (comparison of the 
results calculated for bias voltages in the range Vb = —0.3 
to Vb = —0.1 V show that there are very small vari- 
ations in the length dependence and magnitude of the 
voltage distribution in the system, within the considered 
range of bias voltage). However, as the channel length 
is successively decreased below ss 1.5 /xm, the fraction 

V c hannel 

/V b gradually decreases, while V cont acts/V b grad- 
ually increases, to the point where, for L c < 0.3 /xm, the 
majority of the potential drop occurs across the contact 
regions. Note that the calculated curves for the CC and 
CS size cases are almost identical, showing that no con- 
tact size effects are responsible for the observed features. 

The observed results and dependence of V c hannei/V b 
and V con tactslVb show that a crossover occurs around 
L c w 1.5 /jm and that for shorter channel lengths, bal- 
listic effects become important. With decreasing channel 
length, the transfer of the voltage drop from the channel 
region to the contact regions, indicates that the transport 
becomes increasingly ballistic and it is expected that for 
very short channels, the majority of the voltage drop will 
occur across the contact regions, and only a very small 
(but finite) resistance will occur within the channel re- 
gion. The observed features in the electrostatic response 
are a signature of the crossover between the diffusive and 
quasi-ballistic regimes of transport. 
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FIG. 4: (Color online) (a) Electron densities, n + (thick lines) 
and n~ (thin lines) for the same structure and bias voltages 
as in Fig. [3] Inset shows the minimum densities, ra+ in (solid 
line) and n~ in (dashed line) as a function of VJ,. (b) Minimum 
channel densities, n+ in and n~ in , as a function of channel 
length. Inset shows electron density profiles, n + (solid lines) 
and n~ (dashed lines), as a function of x/L c , for different 
channel lengths, L c . 



2. Electron density 

Next, we perform a similar study of the electron den- 
sities n + and n~ , defined earlier in the text. In Fig.^Ja) 
we show the electron density around the channel region 
for the same structure as in Fig. |^a), calculated at the 
same bias voltages. The thick and thin lines correspond 
to the n + and n~ components of the total density, respec- 
tively. At Vb — the two components are naturally cen- 
tered, although slightly shifted from each other, around 
the channel mid-point. The slight shift is due to the fact 
that, as seen in Fig.[3{a), the electric field is non-zero ev- 
erywhere but at mid-channel, and anti-symmetric with 
respect to the channel center. Hence, the electron distri- 
butions in each half of the channel consist of a slightly 
shifted MB distribution. In the channel half close to the 
source, the distribution is shifted toward negative veloc- 
ities, whereas the distribution in the other half is shifted 
toward positive velocity, see Fig. Eld) for a contour plot 



of the electron distribution around the channel region at 
zero bias. Consequently, ?i + — n~ =/= around the channel 
region (except at the channel mid-point), where the built- 
in field is non-zero. Note, however, that the total electron 
density is entirely symmetric along the x-direction (not 
shown here), as required since 14=0. 

As a bias voltage is applied, two main observations can 
be made. First, the difference between the two compo- 
nents, n + — n~, increases considerably with increasing 
\Vb\, which is also shown in the inset to the figure for the 
minimum densities, TiZ^ and n~- , along the channel. 
Second, the spatial dependence of n + is very different 
from that of n~, and more importantly varies consider- 
ably for different applied voltages. In particular, the po- 
sition of n+ in changes significantly with Vb, whereas the 
position of n~ in is roughly bias independent, as indicated 
by the schematic lines in the figure. This is an indication 
that the efficiency of the scattering is different for the 
v + and v~ electrons. Interestingly, the relative excess of 
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min 



v + electrons in the channel region reaches its maximum 
around mid-channel (not shown), as previously seen in 
Fig. Gib) for Vb — —0.3 V, and the position is relatively 
insensitive to the applied bias voltage. This is another 
indication that the behavior in the channel is controlled 
by the intrinsic fields and doping profiles rather than the 
externally applied electric fields. 

As the channel length is increased (for Vb = —0.3 V) 
the spatial dependence of the electron densities n + and 
n~ changes significantly. In the inset of Fig. Efb) we 
show n + (solid lines) and n~ (dashed lines) for differ- 
ent channel lengths ranging between L c = 0.2 — 4 /Ltm, 
as a function of x/L c . In addition, Fig. HJb) shows the 
dependence of the minimum densities, « min and w min as 
a function of the channel length, calculated at —0.3 V 
for the CS case (we verified again that the contact size 
and length scaling scheme do not influence the results 
qualitatively). 

It is seen that: First, as the channel length increases, 
the electron density in the channel region decreases, ap- 
proaching the value of the background doping, 10 19 
Second, the large magnitude difference between n 
and n~ in , observed for short channel lengths, decreases. 
Third, a crossover from a superexponential to an expo- 
nential dependence of n^^J on the channel length is 
observed around » 1.5 jum. The latter observation we 
again interpret as a signature of a transition between the 
quasi-ballistic and diffusive regimes of transport, and oc- 
curs around similar channel lengths as observed in the 
length dependence of the electrostatic response, discussed 
in Fig.OJb). For long channel lengths, the density depen- 
dence approaches the purely diffusive regime in which 
the n + and n~ have a similar spatial dependence and 
approach the values of the background doping. 

It is interesting to compare the results for the total 
electron densities, obtained by the integration of the elec- 
tron distribution over the entire velocity space, see eq. 
©, with the "ballistic" electron densities, defined in eq. 
. In Fig. we show the fraction of positive and neg- 
ative velocity components, and rig, of defined in 
eq. 0, relative to the total density. 

In Fig. OJa) we show the results for the same 200 nm 
channel length structure studied in Figs. E^a), and 
Ufa), calculated for the bias voltages Vb — —0.1 —0.2 
and —0.3 V. With the decrease of bias voltage, the two 
peaks corresponding to n~^/n (thick lines) and rig/n 
(thin lines) are shifted toward the middle of the chan- 
nel, and decrease and increase, respectively, as expected 
by the smaller asymmetry and magnitude of the electric 
field. We note that, throughout the presented bias volt- 
age range, n~g/n decays approximately linearly on both 
sides of the peak centered at the top of the potential bar- 
rier at the source-drain interface. The fraction of nega- 
tive "ballistic" electrons, rig/n, on the other hand, has 
a faster, exponential-like decay on the drain side of the 
potential barrier, indicating diffusive-like transport. 

The length dependence of the ballistic fraction of elec- 
trons is found to show a distinct signature of the crossover 
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FIG. 5: (Color online) "Ballistic" electron densities, rig and 



n B as denned in the text 



for (a) the same structure and 
applied voltages as in Fig. P|a). and (b) different channel 
lengths, L c , vs. relative channel position, x/L c . 



between the diffusive and quasi-ballistic regimes of trans- 
port, at similar channel lengths where the crossover was 
observed in the electrostatic response as well as the total 
electron density distribution. 

In Fig. Etb), we show n\jn (solid lines) and rig/n 
(dashed lines) as a function of x/L c . Interestingly, n^/n 
has similar decay lengths, relative to L Cl and an almost 
vanishing contribution around mid-channel (with the ex- 
ception of the structure with L c = 0.2 jiva). We attribute 
this behavior to the diffusive nature of the transport of 
electrons with v < 0, since the resistance in the diffusive 
regime is linear with respect to the channel length. 

The positive velocity ballistic electron fraction n~^/n 
on the other hand has: i) a similar decay length and 
length dependence with n^/n for L c — 4 /xm, ii) an 
increasing relative decay length for decreasing channel 
length, yielding a significant contribution throughout the 
full extension of the channel. More importantly, it is in- 
teresting to note the transition between an exponential- 
like dependence for long channel lengths, to an almost lin- 
ear dependence at shorter channel lengths, as a function 
of the relative position, x/L c , within the channel region. 
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We interpret this as yet another signature of the crossover 
between the diffusive and quasi-ballistic regimes of trans- 
port. 



3. Electron distribution function 

Having discussed the signatures of the quasi-ballistic 
and diffusive regimes of transport in the electrostatic 
response and the electron density distribution we now 
direct our attention to the full nonequilibrium distribu- 
tion function. In Figs. El we show the calculated normal- 
ized distribution function, f (x , v)vth / n(x) , for 200 nm 
channel structures for different applied voltages. Fig- 
ure E( a) shows the calculated normalized distribution at 
Vb = —0.3 and is identical to the one shown in Fig.Gfb), 
and thus displays the previously discussed ballistic, high- 
velocity peak throughout the channel region. This peak 
shifts rapidly toward higher velocities along the direction 
of transport and dominates the distribution function well 
into the second half of the channel, close to the drain side 
of the system. In the vicinity of the drain, however, as 
previously discussed, drain-injected electrons spill over 
into the channel and dominate the electron distribution. 
Note however that the two contributions are well sepa- 
rated in phase-space and thus can be easily identified. 

As the bias voltage is decreased, the ballistic signatures 
become less pronounced. Although there still are ballis- 
tic peak-like features for Vb = —0.2 and —0.1 V, see Figs. 
|HJb) and (c), respectively, the maximum velocity and the 
penetration depth of the peak into the channel, i.e., the 
region in which it dominates the electron distribution, are 
decreased. Furthermore, a smoother transition occurs 
between the main peaks in the distribution on the source 
and drain side of the channel [see Fig.HJc)] indicating an 
increasingly diffusive character of the transport. In Fig. 
Eld), for Vb — V, the distribution is as expected sym- 
metric along the direction of transport, but nevertheless 
reflects the inhomogeneous properties of the sample, as 
seen in the slightly shifted distributions around each of 
the two interfaces. 

In Fig. we show the normalized electron distribution 
functions, calculated at Vb — —0.3 V, and for the channel 
lengths L c — 0.2, 0.8, 1.4, and 4 ^m. We again include 
the results for the 200 nm long channel and Vb = —0.3 V 
biased system for comparison in Fig.[7{a) [same as shown 
in Fig. Efb) andEJa)]. Also, for the sake of comparison 
between the different structures, we show the normalized 
electron distribution function as a function of velocity, v, 
and normalized spatial position in the channel, xj L c . 

For L c = 0.8 fim, it is clear that the quasi-ballistic 
features are persistent, although weaker, i.e., lower maxi- 
mum velocity and extension into the channel region (rela- 
tive to the channel length), as well as a smaller separation 
between the drain- and source-injected electron distribu- 
tions. Note again that we plot the results against x/L c 
and that here, the ballistic peak extends over 300 nm into 
the channel. 
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FIG. 6: (Color online) Normalized electron distribution, 
f(x,v)vth/n(x), for the structure with L c = 200 nm and the 
same parameters as studied in Figs. and calculated at 
the bias voltages, (a) V b = -0.3 V, (b) V b = -0.2 V, (c) 
V b = -0.1 V, and (d) V = 0. 

For L c = 1.4 /im, the distribution function shows 
significantly reduced ballistic features. At this channel 
length, there is no clear formation of a ballistic peak 
as in the previous cases. Furthermore, the distribution 
has comparable width, only one dominant peak through- 
out the entire channel region and a smooth connection 
between the near-equilibrium source and drain distri- 
butions. This indicates that the transport is close-to- 
diffusive. For longer channels, the electron distribution 
becomes even more featureless as the transport becomes 
more and more diffusive. In Fig. |7Jd), the distribution 
function for a L c = 4 /jm channel based system is shown 
for comparison. No ballistic features are present, the ve- 
locities are low in comparison to the previous cases and 
the distribution in the channel region, although broad- 
ened and shifted due to the increased fields, is close to 
the electron distributions in the source and drain con- 
tacts. For channel lengths shorter than L C —1A fj,m we 
have, however, seen a transition toward an increasingly 
ballistic regime of transport, with features resembling the 
ones shown in Figs. Ufa) and (b). 

IV. SUMMARY AND CONCLUSIONS 

In summary, we have performed a detailed micro- 
scopic study of quasi-ballistic transport in deep submi- 
cron semiconductor channels. In particular, we studied 
the crossover between the quasi-ballistic and diffusive 
regimes of transport and identified signatures in the elec- 
trostatic response, electron density and full, nonequilib- 
rium electron distribution function, that are due to bal- 
listic transport. The study was performed by numerically 
solving the coupled, nonequilibrium Boltzmann-Poisson 
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FIG. 7: (Color online) Normalized electron distribution, 
f(x,v)vth/n(x), for a structure with externally applied bias 
voltage Vb — —0.3 V, calculated for the channel lengths (a) 
L c = 0.2 pirn, (b) L c = 0.8 /im, (c) L c = 1.4 /im, and (d) 
L c = 4 lira.. 



system of equations, using a self-consistent direct method 
developed by us, that enables us to obtain noise- free elec- 
tron distributions throughout the interesting submicron 
channel regions. 

Semiconductor systems on the deep submicron and 
nanometer scale can in general very easily be brought 
out of equilibrium due to their reduced dimensionality. 
As shown in our model structures, even for moderately 
downscaled channels and realistic doping levels, the built- 
in fields are high and have a very strong spatial de- 
pendence. Electron transport subject to these fields is 
therefore highly out-of-equilibrium and drift-diffusion ap- 
proaches break down. 

Our self-consistent calculations of the electron distri- 
bution function in submicron GaAs channels using the 
Boltzmann-Poisson equations show that the electron dis- 
tribution is radically different from a shifted Maxwell- 
Boltzmann distribution. The distribution is far from 
Gaussian-shaped, has a large broadening, and displays 
pronounced features, peaks and shoulders in the high- 
velocity tail of the distribution. The high-velocity fea- 
tures are signatures of ballistic transport, which occur 
due to the large and strongly spatially varying electric 
fields. The large broadening and significant contribution 
to the negative-velocity component of the distribution 
are, however, in part due to backscattered electrons in 
the channel, and partly due to electrons injected from 
the drain. 

We performed an analysis of the spatial density distri- 
butions for electrons with positive and negative veloci- 
ties, respectively. It was shown that the relative scatter- 
ing efficiency varies within the channel region and gives 



rise to asymmetric density distributions for the forward 
and backward traveling electrons. We also defined and 
studied a "ballistic" electron density, corresponding to 
electrons with total energies higher than the energy of 
the top of the barrier at the source-channel interface. It 
was shown that the spatial dependence for the positive 
and negative velocity contributions to the total "ballis- 
tic" density is radically different. This is due to the large 
difference in the mean free path between the forward 
and backward traveling electrons caused by the strong 
built-in electric field. It was shown that a significant 
number of the source-injected "ballistic" electrons reach 
the channel-drain interface and propagate into the drain. 
The relative fraction of ballistic electrons in the drain is, 
however, very low due to the orders of magnitude higher 
electron density in the drain contact. 

We have shown how these signatures are affected by an 
externally applied bias voltage, and more importantly by 
the length of the channel region. The effects of scattering 
have thereby explicitly been illustrated. 

In general, we have found that a crossover between 
the diffusive and quasi-ballistic regimes of transport for 
the GaAs structures with the given doping concentra- 
tions studied at room temperatures occurs around chan- 
nel lengths L c sa 1.5 /im, which is far longer than the 
mean free path of « 65 nm determined by the ther- 
mal velocity, v t h- The main parameters that influence 
this crossover are the effective mass, the scattering time 
and the details of the doping profile. Hence, a full self- 
consistent treatment such as ours is needed to accurately 
model the crossover between the two regimes of trans- 
port. 

In the electrostatic response we found that the 
crossover manifests as a redistribution of the voltage drop 
between the contacts and the channel region. For large 
channel lengths, the transport is diffusive and the major- 
ity of the potential drop occurs across the channel region. 
For channel lengths smaller than « 1.5 /im a monotonic 
redistribution of the potential drop from the channel re- 
gion to the contacts occurs, and a crossover between the 
two fractions occurs around « 0.3 ^m, beyond which, for 
decreasing channel lengths, the majority of the potential 
drop occurs across the contact regions. 

Similar crossovers were observed in the total and "bal- 
listic" electron density distributions, ti/ + '~) and rig ' \ 
as well as in the full, nonequilibrium distribution, f(x, v). 
In particular, we note that that the spatial distribution 
of relative to x/L c has an increasingly larger decay 
length and more importantly shows a crossover between 
an exponential-like to a linear decay rate around chan- 
nel lengths at which the crossover between the diffusive 
and quasi-ballistic regimes of transport were observed in 
the electrostatic response. The spatial distribution of 
n^, however, has a relatively similar spatial distribution 
for different channel lengths, relative to L Cl and displays 
exponential-like decay indicative of diffusive transport. 
On a microscopic level, we also showed explicitly how 
the signatures of the crossover manifest in the nonequi- 
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librium electron distribution, f(x,v). 

In general, scattering and quasi-ballistic transport in 
ultrasmall semiconductor channels exhibit complex fea- 
tures due to the inherent nonlocal nature of the problem. 
Our study has provided insight into the microscopic sig- 
natures of the problem as well as of the crossover between 
the diffusive and quasi-ballistic regimes of transport. We 
note that some of the features reported in this paper re- 
semble experimental observations that may be explained 
within the presented physical picture. 
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